Proceedings of ASME 2014 Smart Material, Adaptive Structures and Intelligent Systems 

SMASIS2014 

September 8-10, 2014, Newport, Rhode Island, USA 


SMASIS201 4-7645 


GUIDED WAVE DELAMINATION DETECTION AND QUANTIFICATION 
WITH WAVEFIELD DATA ANALYSIS 


Zhenhua Tian 

Department of Mechanical Engineering, 
University of South Carolina, 
Columbia, SC 29208, tianz@email.sc.edu 

Jeffrey P. Seebo 

Analytical Mechanics Associates, Inc. 
Hampton, VA 23681 , ieffrev.p. seebo@nasa.gov 


ABSTRACT 

Unexpected damage can occur in aerospace composites 
due to impact events or material stress during off-nominal 
loading events. In particular, laminated composites are 
susceptible to delamination damage due to weak transverse 
tensile and interlaminar shear strengths. Developments of 
reliable and quantitative techniques to detect delamination 
damage in laminated composites are imperative for safe and 
functional optimally-designed next- generation composite 
structures. In this paper, we investigate guided wave 
interactions with delamination damage and develop 
quantification algorithms by using wave field data analysis. The 
trapped guided waves in the delamination region are observed 
from the wavefield data and further quantitatively interpreted 
by using different wavenumber analysis methods. The 
frequency-wavenumber representation of the wavefield shows 
that new wavenumbers are present and correlate to trapped 
waves in the damage region. These new wavenumbers are used 
to detect and quantify the delamination damage through the 
wavenumber analysis, which can show how the wavenumber 
changes as a function of wave propagation distance. The 
location and spatial duration of the new wavenumbers can be 
identified, providing a useful means not only for detecting the 
presence of delamination damage but also allowing for 
estimation of the delamination size. Our method has been 
applied to detect and quantify real delamination damage with 
complex geometry (grown using a quasi-static indentation 
technique). The detection and quantification results show the 
location, size, and shape of the delamination damage. 
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1. INTRODUCTION 

Aviation safety is an area of focus in NASA’s aeronautics 
work. Damage due to impact events must be detected and the 
damage must be evaluated to determine if repairs are needed 
for safe operation. Matrix cracking, fiber breakage, fiber- 
matrix debonding, and delamination are the most common 
types of damage in fiber composites [1]. In particular, 
laminated composites are susceptible to delamination damage 
due to their weak transverse tensile and interlaminar shear 
strengths (as compared to their in-plane properties) [1]. 
Delamination damage can occur largely internally to the 
composite such that a damage indication is barely visible on the 
composite surface to the naked eye; such damage is known as 
Barely Visible Impact Damage (BVID) [2]. Damage below the 
BVID limit is not visible on the composite surface, even if 
delamination damage exists internally. Such hidden 
delamination damage must be detected and evaluated before it 
becomes critical. The problem of how best to detect such 
material changes in laminated composite components is still 
being solved, and a solution to this problem will be imperative 
for safe and functional, optimally-designed next-generation 
composite aircraft. Current composite “safety margin” design 
comes with a weight penalty which could be mitigated by the 
development of accurate, efficient, and reliable structural health 
monitoring (SHM) methods to locate and quantify the unique 
damage types that can occur in composite materials. 

Guided ultrasonic waves have proven effective for SHM 
due to their ability to inspect a large area while maintaining the 
sensitivity to small defects in the structure [3-8]. Advances in 
guided wave based SHM technologies over the last decade have 
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demonstrated the feasibility of detecting and locating damages 
in composite structural components [9-17]. Due to the complex 
wave interaction that occurs when hidden delamination damage 
is present, tremendous efforts have been put into delamination 
detection methodologies and the study of wave interaction with 
this damage type, such as references [17-23]. 

In this paper, we study guided wave interactions with the 
delamination damage and apply wavefield analysis 
methodologies to analyze the complex wave interaction 
phenomenon and to quantify the delamination damage. A 
scanning laser Doppler vibrometer is used to measure the time- 
space wavefield of guided wave interactions with the 
delamination. Three dimensional (3D) Fourier transform (FT) 
was used to transform the wavefield data from the time-space 
domain into frequency-wavenumber domain where further 
wave propagation information is revealed. To detect and 
quantify the delamination damage, we apply a wavenumber 
field analysis method which can provide the wavenumber value 
at each spatial location. The wavenumber field analysis method 
is then applied to inspect a plate with impact induced 
delamination damage. The detection and quantification results 
show clear indication of the location, size, and shape of the 
delamination damage. 

2. GUIDED WAVES IN LAMINATED COMPOSITES 

When guided waves propagate in a delaminated 
composite, multiple reflections within the delamination region 
can occur, as shown in references [24-28]. Hence, a 
considerable amount of ultrasonic energy is “trapped” above 
and below delaminated regions until eventually the energy 
dissipates. Ramadas et al. studied the interaction of the anti- 
symmetric Ao mode with symmetric delaminations [25]. They 
also identified multiple reflections in the delamination region. 
The trapped energy phenomenon of guided wave propagation 
in delaminated composites has also been studied experimentally 
through the use of wavefield images obtained from laser 
vibrometry tests [26, 27]. It is observed that upon entering the 
delaminated region, the incident waves will be split above and 
below the delamination and then propagate independently 
through upper and lower lamina (single delamination scenario 
being assumed). After that, a significant portion of the waves 
are reflected back from the far edge of the delamination (exit 
region). The reflected waves then pass back above (or below) 
the delaminated area and experience reflections at the original 
entrance to the delamination region. The process is repeated 
until the energy is dissipated, as found in the simulation work 
reported by Hayashi and Kawashima[24]. Glushkov et al. also 
observed similar wave energy trapping and localization in a 
strip delamination by using laser vibrometry to record the 
wavefield [28]. 

Since the material layers above and below the 
delamination have different thicknesses from the original plate 
thickness, the trapped waves will accordingly have different 
wave propagation characteristics. Therefore, it is expected that 
by analyzing the trapped wave signals, the delamination 
dimensions can be quantified. Prior work by Rogge and Leckey 


demonstrated that wavenumber analysis can be used to 
determine the approximate depth and size of near-surface 
delamination damage [23]. This paper demonstrates 
quantitative characterization of the trapped waves, along with 
delamination detection and sizing utilizing the trapped wave 
phenomenon for real delamination damage. 

To demonstrate the fundamentals of guided wave 
propagation and interactions in laminated composites, two 380 
mmx380 mm 8-ply IM7/8552 composite laminate plates with 
layup [02/902 ]s are first studied. One plate was pristine, while 
the other contained a delamination created by inserting a 
20mmx20mm Teflon film under the second ply during 
manufacturing. Guided waves are actuated by a surface 
bonded 7 mm round PZT patch at 300 kHz, while the time- 
space wavefield data are acquired by a scanning laser Doppler 
vibrometer and are used to investigate the wave propagation 
[29]. Figure 1 shows the schematic of the test plate (detailed 
experimental setup can be found in [29]). 

Figure 2 shows the time-space wavefields of the guided 
waves propagating along the y-axis (ID wave propagation). 
For the pristine case in Figure 2a, a faster So mode and a slower 
Ao mode can clearly be observed. For the delaminated case in 
Figure 2b, waves propagating back and forth between the 
beginning (y=20 mm) and ending (y=40 mm) boundaries of the 
delamination can be distinctively seen, exhibiting a “trapped 
wave” phenomenon [27]. 

Figure 3 shows the snapshots of the wavefields in the scan 
area. Figure 3 a is the wave propagation recorded at 15 ps when 
the incident So mode interacts with the delamination damage. 
We see there are waves propagating both upward and 
downward between delamination boundaries (upper and lower). 
At 30 ps, as shown in Figure 3b, the incident Ao mode is 
interacting with the delamination damage. The figure shows a 
localized yet complicated interference pattern in the 
delamination region due to the multiple trapped waves and 
wave interferences at this time. 



PZT actuator 


Unit: mm 


Figure 1. Composite piate setup for the iaser vibrometry 
test. 
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Figure 2. The time-space wavefieids of the waves 
propagating aiong the y axis: (a) for the pristine piate, and 
(b) for the deiaminated piate. The two soiid biack iines at 
20mm and 40mm are the deiamination boundaries. 



(a) X coordinate (mm) 



Figure 3. Guided wave propagation and “trapping” in the 
deiaminated composite piate: (a) the So mode interacting 
with the deiamination at 15 ps, and (b) the Ao mode 
interacting with the deiamination at 30 ps. The biack soiid 
iines are the deiamination boundaries. 


3. DATA ANALYSIS APPROACHES 

Though interesting wave interactions have been observed 
in the time -space wavefield, a quantitative evaluation is not yet 
readily obtained. To quantify the complicated wave interactions 
within the deiamination region, further wave data analysis 
becomes essential. Wave propagation characteristics such as 
wave mode content, as well as how the waves evolve as with 
frequency or space are also desired. Prior works reported in the 
literatures have demonstrated that wavefield data in the 
frequency-wavenumber representation have abundant 
information regarding the existence of various wave modes and 
wave propagation characteristics [21, 30-33]. We therefore 
develop and apply several quantitative wave analysis 
approaches based on wavenumber information in order to study 
wave interactions with deiamination type defects, and to detect 
deiamination damage quantitatively. 

3. 1. Frequency-wavenumber spectral analysis 

Frequency- wavenumber information of the guided waves 
can be obtained from the time-space wavefield by applying a 
3D FT, mathematically given as 


= j I u{t,x)e (1) 

where the space vector x and the wavenumber vector k in 
Cartesian coordinates are defined as (x, y) and (fc, ky). u{t,x) is 
the time-space wavefield in terms of the time t and the space 
vector X. U{f, k) is the resulting frequency-wavenumber 
representation or “spectrum” in terms of the frequency / and the 
wavenumber vector k. As the frequency is the counter part of 
time, the wavenumber is the counter part of space. The 
mathematics of the frequency-wavenumber technique can be 
found in several references [21, 30-34]. 

Figure 4a and 4b show plots of the frequency- 
wavenumber spectra of the experimental wavefieids for the 
pristine and delaminated cases, respectively. To identify the 
wave mode information, the frequency-wavenumber spectra are 
compared with the theoretical dispersion curves which are 
plotted as dotted (Ao) and solid (So) lines. The pristine plate 
frequency- wavenumber spectrum (Figure 4a) shows the lower 
amplitude Ao and higher amplitude So modes, which match well 
with the theoretical dispersion curves. 
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In the delaminated plate, the delamination damage splits 
the originally bonded full laminate into two parts (above and 
below the delamination). Since the two separated parts have 
different layups and thinner thicknesses compared to the 
original pristine laminate, the guided waves in the separated 
parts will have modified propagation characteristics (such as 
modified dispersions and wavenumbers). 

As expected, the spectrum of the delaminated case (Figure 
4b) shows new frequency-wavenumber components in addition 
to the original Ao and So modes. The new frequency- 
wavenumber components include both positive and negative 
counterparts, which correspond to the forward and backward 
propagating waves, respectively. We compared the new 
frequency-wavenumber components with the dispersion curves 
of the material above the delamination, equivalent to a two ply 
[0\0] laminate, since only the waves in the top plate surface was 
measured in the experiment. The comparison in Figure 4b 
shows the new frequency- wavenumber components match well 
with the theoretical curves of Ao mode in the two ply [0\0] 
laminate. This analysis result agrees well with the real 
delamination depth which is below the second ply. 


3.2. Spatial wave filtering 

The frequency- wavenumber analysis result shows that 
new frequency-wavenumber components appear in the 
spectrum of the delaminated plate, and are correlated to the 
waves in the delamination region. To further confirm this 
result, a wave filtering process is used to extract the new 
frequency-wavenumber components and transform them back 
into the time-space domain. The filtering process is 
mathematically expressed as the product between the 
frequency-wavenumber spectrum U(f, k) and a filter function 
F(f, k) as: 

t/^(/,k) = t/(/,k)F(/,k) (2) 

where F(f, k) denotes a filter in the frequency-wavenumber 
domain and Ufif, k) is the filtered spectrum, which only keeps 
the delamination induced new wavenumbers. By transforming 
the filtered spectrum Upif, k) back into the time -space domain 
through inverse 3D FT, we can obtain the filtered time-space 
wavefield Ufit, x) as 

= (3) 


(a) 


(D 

§ 



Frequency (kHz) 



Figure 4. The frequency-wavenumber spectra of the time- 
space wavefieids in Figure 2: (a) for the pristine piate, and 
(b) for the deiaminated piate. The soiid and dotted iines are 
theoreticai dispersion curves for So and Ao wave modes, 
respectiveiy. 


Since the filtered spectrum k) only keeps the new 

frequency- wavenumber components, it is expected that the 
corresponding time-space wavefield up{t, x) only contains the 
trapped waves induced by the delamination. 

Figure 5 shows the filtered time-space wavefield along the 
y axis. Comparing this result to the wavefield before filtering 
(Figure 2b), it is seen that the filtered time-space wavefield 
shows trapped waves that exist between the delamination 
boundaries. 

Figure 6 shows a snapshot at 3 Ops of the filtered 
wavefield in the scan area. The filtered wavefield only contains 
the trapped waves in the delamination region, and shows the 
feasibility of using the trapped waves for delamination size 
quantification. Moreover, the wavefield filtering results further 
confirm that the new frequency-wavenumber components in 
Figure 4b correspond to the trapped waves in the delamination 
region. 



0 10 20 30 


Time (ps) 


Figure 5. The filtered time-space wavefield of the waves 
propagating along y axis. The two solid black lines at 20mm 
and 40mm are the delamination boundaries. 
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X coordinate (mm) 


Figure 6. A snapshot at 30ps of the filtered wavefield in 
the 2D scan area. The black solid lines are the delamination 
boundaries. 


3.3. Wavenumber field analysis 

The frequency-wavenumber analysis presented in the 
section 3.1 unveils wave propagation characteristics that cannot 
be explicitly seen in the time-space domain. However, in the 
resulting frequency-wavenumber representation U(f, k), the 
spatial information x is not available from the result of the 3D 
FT (see Figure 4). 

It is quite often desirable to know how the wavenumber 
changes in space. To retain the spatial information, a novel 
short space 3D FT has been applied to acquire a space- 
frequency-wavenumber representation. The technique is a 
straightforward extension of the short time FT to multi- 
dimensional problems, i.e., breaking down the time-space 
wavefield into small segments over the space dimension before 
Fourier transformation [35]. To do this, the wavefield data is 
multiplied by a window function which is non-zero for only a 
short period in space while constant over the entire time 
dimension. The 3D FT is then applied to the resulting wavefield 
segments. When the window slides along the space dimension, 
a set of windowed wavefield segments is generated. The 3D FT 
is applied to these segments, resulting in a set of frequency- 
wavenumber spectra that are indexed by the locations of the 
window. Through this technique the spatial information is 
retained. Mathematically, the short space 3D FT is 
implemented through spatial windowing technique, as 

S(X,/,k) = f" f” (4) 

J ^X) J -CO 

where X is the retained spatial vector and W(t, x) is a 
window function. In our study, a Hanning function is used to 
construct the window W(t, x), given as: 




( Ixl Y 



0.5 

1+cos 

A. 

In — 
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VI 

(5) 


0 otherwise 



X coordinate (mm) 


Figure 7. The calculated wavenumber field from the 
experimental time-space wavefield in the delaminated plate. 
The black solid lines are the delamination boundaries. 


where Dx is the window length in the space domain. Further 
details of the short space 3D FT can be found in [23, 32, 33]. 

As stated above, by sliding the window W{t, x) along the 
space dimension, the resulting space-frequency-wavenumber 
spectrum 5'(X,/,k) can indicate how the frequency and/or 
wavenumber components vary in space. Thus, by finding the 
wavenumber vector k* with the highest spectrum amplitude 
for each space and frequency point, we obtain the wavenumber 
vector function k'"(X,/) , as 

k’(X,/) = argrmx|5'(X,/,k)| (6) 

Then, for a specific frequency fo , we define the wavenumber 
field as 

^•(X) = |k*(X,/„)| (7) 

For a series of frequencies fi (/=1,2,3...A) within a small 
frequency band, the wavenumber field is defined as 

^’(x)=T|]|k*(Xy;)| (8) 

The wavenumber field A:*(X) contains the local wavenumber 
information at each spatial point. Figure 7 gives the 
wavenumber field calculated from the time-space wavefield of 
the delaminated plate. The wavenumber field shows an area 
which has larger wavenumber values than the rest. This region 
matches well with the real delamination location and size. 
Since the delamination damage affects the wavenumber (as 
demonstrated in section 3.1), the wavenumber field can be used 
to detect and quantify the delamination location and size. 
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4. IMPACT-INDUCED DELAMINATION DETECTION 
AND QUANTIFICATION 

The previous section has shown that the wavenumber 
analysis approaches yields quantitative information of guided 
wave interactions and wave trapping in delamination regions; 
and can also be used as delamination detection tools. In this 
section, the wavenumber analysis methods are used to detect 
and quantify real impact-induced delamination damage in 
laminated composite plates. 

4.1. Damage growth in a CFRP plate 

A quasi-static indentation technique was used to grow 
damage in a CFRP plate. The technique has been shown by 
previous authors to be an ideal method for controlled growth of 
impact-like damage in composite laminates [3 6]. The composite 
sample is a 26 ply IM7/8552 381 mm x 381 mm laminate with 
a quasi-isotropic layup of [(02/452/-452)2 90] s. The composite 
plate was loaded in a mechanical testing system using a 50.8 
mm diameter indenter at a constant displacement rate of 1.27 
mm/min. The displacement and load were monitored during the 
indentation test and the composite was unloaded when the data 
showed a sudden drop in load (due to damage growth), as 
shown in Figure 8b. A photo of the CFRP plate after the 
indentation test is shown in Figure 9a. Figure 9b is a zoomed in 
view which shows a very slight indentation on the surface with 
no other surface damages. 


(a) 




Figure 8. Quasi-static indentation for deiamination 
damage growth: (a) the indentation setup, and (b) the curve 
of indentation force versus dispiacement. 





Figure 9. The 26 piy iaminate after the quasi-static 
indentation: (a) a photo of the specimen, (b) a zoomed view 
of the indentation, and (c) a Cscan image showing the 
interior deiamination damage generated by the indentation 
test. 


Following the quasi-static indentation test, an ultrasonic 
Cscan was used to quantify the resulting deiamination damage 
under the plate surface. The Cscan data were collected in an 
immersion tank using a 10 MHz transducer and a scan step size 
of 0.25 mm. The Cscan image in Figure 9c shows multiple 
delaminations within the 32.3 mm x 22.6 mm damage zone, 
occurring at various depths through the plate thickness. 

4.2. Deiamination detection and quantification 

A 1 mm round PZT was installed on the CFRP plate 
surface (shown in Figure 9a) to generate guided waves at 500 
kHz. A scanning laser Doppler vibrometer was used to measure 
the wavefield in the CFRP specimen containing the impact- 
induced deiamination damage. Figure 10a is a wavefield 
snapshot at 55 ps showing the guided wave propagations and 
interactions. The incident guided waves interact with the 
damage and are trapped in the deiamination region. Due to the 
existence of multiple wave modes and complex wave 
interactions, it is difficult to detect and quantify the 
deiamination damage only from the original measured time- 
domain wavefield. Further wavefield processing is required. 

To detect and quantify the damage, both the wavefield 
filtering and wavenumber field analysis methods were used. 
Figure 10b gives the wavefield filtering result which only 
retains the trapped waves. The trapped waves clearly indicate 
the deiamination location. 

Figure 10c shows the wavenumber field over the test area. 
The deiamination damage is clearly observed as larger 
wavenumbers in wavenumber field. The detected deiamination 
size in the wavenumber field is approximately 30 mm x 20 
mm, which agrees well with the Cscan result. Moreover, the 
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delamination shape in the wavenumber field closely matches 
the delamination shape in the Cscan image Figure 9c. 
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Figure 10. The delamination damage detection and quantification results: (a) the wavefield snapshot at 55 ps shows the 
guided wave interactions and trappings in the impact-induced delamination damage, (b) the filtered wavefield snapshot at 55 
ps shows the trapped waves at the damage, and (c) wavenumber field shows the location and shape of the damage. 


5. CONCLUSIONS 

This paper presents our studies on guided wave 
interactions with impact induced delamination damage and 
damage quantification techniques for laminated composite 
plates. The wavefield data analysis confirms the wave trapping 
phenomenon. Quantitative interpretation of the wave behaviors 
in the delamination region via different wavenumber analysis 
methods is demonstrated. Among the methods studied, the 
frequency-wavenumber analysis shows that the frequency- 
wavenumber spectrum can indicate the presence of new 
wavenumbers. The reconstructed wavefield images created 
from the new wavenumber components demonstrate that the 
new wavenumbers correlate to trapped waves in the 
delamination region. Additionally, the wavenumber field shows 
how the wavenumber changes as a function of the propagation 
distance. Thus, the location and spatial duration of the new 
wavenumbers can be identified, providing a useful means not 
only for detecting the presence of delamination damage but also 
sizing the damage. Finally, the wavenumber field analysis 
shows the wavenumber values at each spatial point indicating 
not only the location and size, but also the shape of the damage. 

The developed analysis methods have been applied to 
detect impact-induced delamination damage in a carbon fiber 
reinforced polymer plate. The delamination is barely visible 
even from careful visual inspection while C-scan images show 
internal delamination of about 30 mm x 20 mm size. By using 
wavenumber field analysis, interior delamination damage is 
detected and quantified. 
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